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Abstract 

We study a family of Hamiltonians of fermions hopping on a set of lattices in the presence of a background gauge 
field. The lattices are constructed by decorating the root lattices of various Lie algebras with their minuscule rep- 
resentations. The Hamiltonians are, in momentum space, themselves elements of the Lie algebras in these same 
representations. We describe various interesting aspects of the spectra — which exhibit a family resemblance to the 
Dirac spectrum, and in many cases are able to relate them to known facts about the relevant Lie algebras. Interest- 
ingly, various realizable lattices such as the kagome and pyrochlore can be given this Lie algebraic interpretation and 
the particular flux Hamiltonians arise as mean-field Hamiltonians for spin-1/2 Heisenberg models on these lattices. 



1. Introduction and Outline 

In this paper we study a family of Hamiltonians of fermions hopping around on various lattices in the 
presence of a background gauge field. These Hamiltonians are interesting to us, and we hope to the reader 
as well, in three distinct contexts: the search for flux phases in quantum magnets, the theory of Lie algebras, 
and their possessing interesting continuum limits with a family resemblance to the Dirac Hamiltonian. Let 
us now expand on these connections. 



1.1. Flux phases 



The first of these contexts is where we encountered them, which is the search for flux phases in quantum 
magnets. It was first noted by Baskaran, Zou, and Anderson j 1 1 12] that a novel mean-field theory for SU(2) 
invariant spin-1/2 Hamiltonians could be constructed by re-representing spins as fermionic bilinears, S = 
ipcnp, and relaxing the constraint of a unit fermionic occupation of each site in favor of a global constraint 
of a half filled band. The mean-field treatment consists of replacing the starting Hamiltonian, quartic in 
the fermions, by one quadratic in them wherein the fermions hop on the lattice in the presence of a self- 
consistently calculated background (frozen) gauge field. Following this, Affleck and Marston [5] noted that 
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by extending the model to SU(N) spins this mean-field theory could be made exact at N — oo whence 
it could serve as the starting point of a 1/N expansion. Specifically, they showed that the Heisenberg 
model on the square lattice exhibited an enticing mean- field solution with a flux of tt (a factor of e Z7r ) per 
plaquette. Remarkably, they found that the solution to this hopping problem led to a Dirac fermion at zero 
energy. Since the system had particle-hole symmetry, and was at half-filling due to its insulating parentage 
mentioned above, the Dirac fermion was right at the Fermi energy and thus central to low energy physics. 
In particular, fluctuations around the saddle point were described by the Dirac fermion minimally coupled 
to the fluctuating gauge field. While it is tangential to our purposes in this paper, we note that the latter 
problem has been the focus of much progress in recent years [4|. Following these early developments, there 
has been much work examining various "flux phase" mean-field theories on various lattices which is too large 
a literature for us to review here0 Most interesting for our purposes was the generalization to time reversal 
(T) breaking phases made by Wen, Wilczek and Zee [9] which thus gave a mean-field meaning to the chiral 
flux phases proposed previously by Kalmeyer and Laughlin [10] on the basis of an inspired ansatz. These 
phases exhibit fluxes through closed loops which are different from the two T-invariant values and w. There 
is one last aspect of this body of work that is also worth noting at the outset, namely that it does not take 
the actual energetics at N — oo, and thus the relative stability of various mean-field solutions, too seriously. 
Strictly at N — oo the kinds of solutions discussed above lose out to fully dimerized mean-field states of 
lower energy [llll2j . Their continuing interest has to do with the possibility that this relative ordering of 
energies is reversed as N is decreased. Indeed, the relative ordering between non-dimerized solutions could 
also change as N is decreased and hence one can (and people do), in good conscience, start out by studying 
various interesting mean-field solutions which at least exhibit local stability [5]. 

With the above recital we can now locate our Hamiltonians: they arise, with a few caveats and excep- 
tions, as mean-field Hamiltonians for nearest-neighbor spin-1/2 Heisenberg Hamiltonians on the appropriate 
lattices. They also involve, in almost all cases, T-breaking. 

1.2. Lie Algebras and Hamiltonians 

There is, of course, a large set of such Hamiltonians and we next need to describe the restrictions that 
generate the family that we study. This brings us to the second context in which our Hamiltonians can 
be situated and which intrigues us most: the deep ties our Hamiltonians have to Lie algebras. These ties 
are twofold: ideas from Lie algebras are central to generating the very lattices the fermions move on and 
the Hamiltonians can be written as direct sums of pieces that are Lie algebra elements in specified repre- 
sentations. Consequently we find that properties of Lie algebras also control some striking features of the 
resulting spectra. Sometimes we fully understand these connections and sometimes we do not, though in all 
cases we will share what we know with the reader. 

Let us first summarize how the unit cell and the underlying lattice have direct group theoretic significance. 
A discussion of the group theory used here can be found in |13|14j . 

Consider the unit cell. Recall that the generators of a (semi-simple) Lie algebras can be partitioned into a 
maximally commuting Cartan subalgebra Hi : [i = 1, . . . r] = H whose eigenvalues label the weights, and a 
set of ladder operators E& and their adjoints E^ = E-a. that act on the states to raise (lower) the weights 
by a. The vectors a are called the roots. The states within any irreducible representation (multiplet) of a 
Lie algebra may therefore be visualized as a collection of points in a space of dimension r, called the rank. 
The coordinates of the points are the simultaneous eigenvalues of H. The roots that help us move around 
these points are also vectors in the same space. For example, in the case of the rank-2 group SU(3), (whose 
commuting quantum numbers are traditionally called isospin and hypercharge in the physics literature) 
the fundamental (quark) representation is an inverted triangle, the anti-quark is a triangle, and the eight- 
dimensional adjoint representation is a hexagon with two null weights at the center. The six nonzero roots 
correspond to the six corners of the hexagon. 



1 See, for example, \E\, [6], [7] and references therein for more recent work in 2 d. A 3 d example and connections to the 
quantum Hall effect were addressed by [§]. 



2 



The unit cells of the lattices we consider correspond to certain special representations called minuscule 
representations. Starting with any one state in a minuscule representation or multiplet, we can obtain all 
the others by acting with the Weyl group, the group of reflections about the hyper-planes normal to the 
roots. All weights of a minuscule representation are on the same footing and in particular all have the same 
length. Thus the quark and anti-quark are minuscule while the adjoint representation, with weights of both 
zero and non-zero length, is not. 

This unit cell is now used to decorate a lattice, which is a subset of the root lattice Lr. Recall that roots, 
like weights, also live in r dimensions. So it is possible to choose r roots, called simple roots, as a basis. A 
basic result of group theory is that every root is an integer linear combination of the simple roots. There are 
of course a finite number of them in any algebra of rank r. The root lattice Lr is the infinite lattice formed 
with the same basis with but with any integer set of coefficients. For SU(3), the six roots form a hexagon, 
while the root lattice is the infinite hexagonal lattice. 

The lattice on which our fermions move is L<ir, the subset of Lr whose points have even integer coefficients, 
decorated by a basis corresponding to a minuscule representation. 



When applied to the quark representation of SU{3), this yields the kagome lattice, as shown in Fig. 1(b) 



The inverted triangle (Fig. 1 (a) ) is the fundamental quark representation which forms the unit cell that 
decorates the root lattice L^r, which is a hexagonal lattice with twice the lattice spacing as the root lattice 
Lr. The origin of coordinates is at the point named q in the figure. If you stare at the figure hard enough, 
you can also see it as a hexagonal lattice decorated by the conjugate representation, the anti-quark, whose 
weights are the negatives of the quark. The center of one such unit cell is labeled q in the figure. The quark 
and anti-quark unit cells are corner sharing. These features are common to all our models and exist because 
the unit cell and the lattice are constructed from weights and roots in a particular way. A proof of this 
will be given in Section 12.11 Let us now turn to the construction of the Hamiltonians on the above set of 
lattices. Gauge fields will enter our models in the form of purely imaginary hopping amplitudes which can 
be H. This restriction means that on any triangular face the flux can only be ±tt/2o In other words the 
background gauge field is an Ising-like variable, and time reversal symmetry is broken. Also, we will require 
that the gauge fields exhibit the periodicity of the Bravais lattice L2R — this has the gauge invariant content 
that there is no net flux passing through the lattice. It is worth noting here that generically in problems 
of this kind we must view symmetries as projective, i.e., the underlying group operations will have to be 
accompanied by additional gauge transformations to make the symmetry manifest [5], say the way Lorentz 
transformations have to be accompanied by gauge transformations in relativistic field theories to establish 
Lorentz covariance. The classification scheme relevant to our problem appears to be that of Color Groups, 
in which each face of the crystal is colored black or white, which we may read as ±7r/2 of flux [15116] . 

We are now in a position to specify the Hamiltonians of interest. Since the flux added to this lattice will be 
translationally invariant on L2R, we may go to momentum space to solve for the dispersion relation. Evidently 
H (k) will be a matrix that acts on the states of the minuscule representation since they constitute the unit 
cell. With our choices of background gauge fields we then arrive at manifestly hermitian Hamiltonians of 
the form 

H(k)= ]T C a (k)(E a + El) (1) 
where the coefficients Cot are real, satisfy 

Ca(-k) = -C a (k) (2) 

and S + are the positive roots. The roots may be divided into positive (£+) and negative (£_) roots by 
drawing a plane through the origin which does not contain any roots. This choice is basis-dependent. For 
any choice of basis we may choose an ordering of the basis vectors so that the positive roots are those whose 
first nonzero component is positive. 

Since the ladder operators move us around the multiplet, it is reasonable to consider H of the form (|TJ) . 
However we must bear in mind that this is not the most general possibility on this lattice. For example the 
model only allows hops between sites that differ by a single root while there are minuscule representations 



Evidently, the flux is the gauge invariant variable. Our choice of gauge fields is convenient for the purposes of this paper. 
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(a) 



(b) 



Fig. 1. The kagome lattice and SU(3): The unit cell (a) is the quark triplet, with weights as indicated. The kagome 

lattice (b) is formed by decorating L2R (the subset of the SU(3) root lattice with even integer coefficients) with this unit cell. 
L2R is a hexagonal lattice with basis vectors ai and &2 which are twice the simple roots of SU(3). The origin of coordinates 
is marked q. The lattice may also be viewed as being decorated by anti-quark unit cells, one of which is centered at q 

where the states differ by more than a single root. Other restrictions implied by this form of H will be 
discussed later. 

Observe that we have obtained an unusual connection between the lattice and the hopping Hamiltonian: 
the former is obtained by decorating the root lattice of a Lie algebra with the weights of one of its repre- 
sentations, and the hopping Hamiltonian (in momentum space) is an element of the same Lie algebra in the 
same representation! It is worth emphasizing that this connection does not involve the symmetry group of 
the problem. For example, the Lie algebra SU(3) that shows up on the kagome lattice does not generate 
the actual symmetries of the hopping problem or even correspond to any symmetries of the starting SU(2) 
Heisenberg model on that lattice. 

One final point on the construction of our lattices and Hamiltonians: Although our lattice construction 
works in all dimensions, we shall limit ourselves to lattices in d = 2 and d — 3, both because these are the 
dimensions we could encounter in the lab and because the rules for attaching flux break down in higher 
dimensions wherein areas cannot be oriented unambiguously. In this latter aspect the flux Hamiltonians 
that give rise to Dirac problems are special as they involve n flux, which does not carry an orientation, and 
thus their generalization to arbitrary dimensions is straightforward. 

1.3. Dirac-like continuum limits 

The discretized Laplace operator on lattices typically looks like a hopping problem with zero background 
flux. It is an interesting fact that when Kogut and Susskind [17] set out to discretize the Dirac operator 
on cubic lattices, they were naturally led to introduce a background flux of tt per plaquette. Indeed, the 
Affleck-Marston work ended up "rediscovering" this earlier construction in the case of two spatial dimensions. 

The flux Hamiltonians we consider evidently have a family resemblance to the 7r flux Hamiltonians which 
they generalize. It is also the case that their low energy limits exhibit a family resemblance to the Dirac 
theory which is the third context in which they appear to be interesting. 

As a first step in explaining what we mean let us observe that zero energy plays a special role in our entire 
set of Hamiltonians. Normally, this comes about via a particle-hole symmetry, where at every momentum k 
states at energy E are accompanied by states at —E. In our problems, the choice of background flux ensures 
that H(— k) = — -ff(k). Consequently for every level that is negative and hence occupied at k, there is one 
at — k that is empty and unoccupied so that the combination of the two bands is particle-hole symmetric 
and E = is again special. We should observe that this is very useful in the original context of the mean- 
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field theory of various Heisenberg models for this ensures that the half-filled band exhibits a Fermi energy, 
E F = 0. 

In the 7r flux phase, the structure about E — is that of the Dirac theory. In our examples we find 
generalizations that we term Dirac-like or pseudo-Dirac. By Dirac-like, we mean a generalization in three 
respects. First, although the spectrum is linear in momentum for small momentum, it possesses only discrete 
and not full rotational invariance. Second, the square of the Hamiltonian H is not proportional to the unit 
matrix, though a higher degree polynomial is. While a result of this sort is inevitable for any finite size 
matrix, the fact that the polynomial often contains just even powers leads to a result that is sufficiently 
reminiscent of the Dirac case. The third generalization we encounter is that in addition to isolated Dirac (or 
Dirac-like) points, wc often find entire lines and even planes of zero energy. 

The zeros arc interesting in and of themselves in the Lie algebraic setting: in some cases, the lines of zeros 
are along the direction of the weights while in other cases the planes of zeros are simply related to the roots 
and so on. In some cases we can understand the locus of zeros without explicit computation by appealing 
to ideas from group theory, while in many cases we could neither anticipate nor explain the zeros. 

There is, inevitably, a matrix structure that goes with the low energy Dirac-like theory and generalizes 
the gamma matrices but it does not appear to be immediately interesting in and of itself although we will 
exhibit it in one especially interesting case. 

1.4. What follows 

In Section [2] we explain the construction of the lattice in detail, and describe the examples we will study 
in the remainder of the paper. We discuss the hopping problem for the d = 2 lattices in Section [3J studying 
several possible background fluxes. In Section|4]we carry out the same analysis for lattices in d = 3. Here the 
reader can find in some detail a discussion of the pyrochlore lattice, which we understand the best in terms 
of group theory and which also displays interesting mathematical structures. In Section [5] we comment on 
the utility of our Hamiltonians as mean-field solutions of the Heisenberg model on our various lattices and 
on two extensions of our analysis. Section [6] contains concluding remarks and discusses open problems. 

2. Lattice construction 

Recall our recipe for generating the lattices: 

- Generate L2R, the even sector of the root lattice, that is to say, even integer combinations of the simple 
roots. 

- Decorate each lattice point with a minuscule representation. 

As noted in the Introduction, although our lattice construction works in all dimensions, we shall limit 
ourselves to lattices in d = 2 and d = 3, due to their physical relevance and more absolutely because the 
rules for attaching flux break down in higher dimensions wherein areas cannot be oriented unambiguously. 
Hence we study rank 2 and 3 groups only. Here is the list of candidates. 

- d = 2; SU(3) and SO(5)=Sp(4). The exceptional group G2 does not have minuscule representations. The 
group SO(4) factors into two independent SU(2) factors and will only be discussed very briefly. 

- d= 3: SU(i)=SO(6), SO(7), and Sp(6). There are no exceptional groups of rank 3. 
The minuscule representations in each case will be listed as we go along. 

2.1. Properties of our lattices. 

The lattices we manufacture by the rules listed above have some interesting features that will be established 
in this section. Before we do so in general, let us pause to examine a simple example, the rank-2 group SU (3). 
This exercise will help us better motivate and understand the general case. 

The only minuscule representations are the quark and antiquark. Let us begin with the quark represen- 
tation. 
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The weights, numbered 1, 2,3 in Fig. 1(a) are 



fx = M 2 =^(-l,^), M, = (0,-^) (3) 



and point to the vertices of an equilateral triangle. 

One significant feature of these weights that we will invoke is that 




1 yjb v — / 

Mi • M< H t . /. (4) 



where a = 1/3 and 6 = —1/6. In other words, there are just two possible values for \i i ■ /Xj, between a weight 
and itself (a) and between a weight and any other (b). 

The nice thing about SU (3) is that in the quark representation, the weights form a simplex where every 
corner is equidistant from every other, the difference between any two corners (any side of the simplex) is a 
root, and these are the only roots. Note that this implies six roots for SU(3), since each edge of the triangle 
can be traversed in two directions. So the root system of SU (3) is 

£(St/(3)) = ^ - = a« [i^j: 1,2 or 3] (5) 

We may choose the simple roots to be a 12 and a.23 and in terms of them the lattice L2R is defined as the 
set of points 

2R = 2mai2 + 2na.ni 

3 (6) 

= msLi + rta 2 



We can now put the expanded root lattice together with the minuscule representation: in Fig. 1(b) we 
have placed one triangle at the origin (labeled q) and made copies at every lattice point in L%r. As the 
reader can see, we end up with the kagome lattice. 

Let us note that the two features, Eqs. (|4l5p . generalize in the obvious manner for all SU(N), where the 
weights now point to the vertices of an TV— simplex. 

Now we turn to some general features of our lattices valid for all the groups we will study, not dependent 
on the special features of SU (N) alluded to above. To see what they might be, look at Fig. l(b)| and observe 
the following features: 

(i) The lattice can also be viewed as L2R decorated by the conjugate (anti-quark) representation with 
reversed weights. The original representation and the conjugate share corners, and every site is shared 
in this manner. 

(ii) If the particle can hop to any point labeled i from a point labeled j in the same unit cell by moving a 
displacement dy , it can keep moving an extra dy to reach a point labeled j in the adjacent unit cell. 
In other words the edges of the unit cell and the conjugate unit cell that meet at a shared corner are 
continuations of each other with no change in direction. 

We will now furnish the proofs of these results in the general case. 

Theorem I: The original lattice 2R + fi i can be rewritten as 2R + 2/i^ — fi t . In other words, if at a new 
origin displaced from the old one by 2fi 1 , we place the conjugate representation and make copies of it using 
any element of L2R, we get the old lattice. The result is just as valid if we use any other weight 2fij in place 
of 2fi 1 . 

Proof: 

2R + 2ji x ^ = 2R + 2/Lt! - 2fi t + i i l = 2R' + fi. (7) 

where we have used the fact that 2/jl 1 — 2fi i; being an even integer multiple of weight differences, is then an 
even integer multiple of roots, which in turn is a translation within i 2 fl- 

Note that the choice of origin at 2/j, 1 is arbitrary: the choice 2/x 2 differs by 2/x 2 — 2^jl 1: an even integer 
multiple of roots, and hence a translation within L,2R- I 

Theorem II If the particle can hop to any point labeled i from a point labeled j in the same unit cell by 
moving a distance d^-, it can keep moving an extra to reach a point labeled j in the adjacent unit cell. 
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Proof: Since dy is a difference of weights it is some integer combination of simple roots. Moving an extra 
distance dy, corresponds to a total displacement by an even integer combination of simple roots, which is 
a symmetry of L2R. It follows that if we start at a point labeled j we must end up a point also labeled j. I 

2.2. Lattices in d = 2 

We have already discussed SU(3) in the last section. Now we will deal with 50(5) and Sp(4). These two 
Lie algebras are mathematically equivalent up to cosmetic differences which will be displayed. 

2.2.1. SO(5) 

We begin with the more familiar group SO(5) which preserves the norm 



i=l 



(8) 



and has a defining representation of 5 x 5 orthogonal matrices. 

We choose as Cartan generators Hi = L\i and Hi = L34 which generate rotations in the 12 and 34 planes 
In terms of the coordinates 

± X\ ± ix 2 j. £3 ± ia^ 



*z - ^ 11 
we may write the invariant in this spherical basis as 



V2 



Xq = X5 



II 



2 ^ x "3jq 



(9) 



(10) 



The vector x itself serves as a 5-dimensional representation. The components xf and are eigenstates 
of Hi, H2 with eigenvalues H = (±1,0) and (0,±1) while xq = x$ does not respond to either rotation and 
has eigenvalues (0,0). The vector representation is not minuscule since the weights are of unequal length. 

The only minuscule representation is the 4-component spinor, with weights 



1 1 
±-,±- 

2' 2 



(11) 



These form a square as in Fig. 2(a) 




Fig. 2. Planar pyrochlore and 50(5): The unit cell (a) of the spinor representation is the square of edge unity. The resulting 
square lattice, with hoppings along all roots, is shown in (b). The lattice vectors ai, a2 are a basis for the subset of the SO(5) 
root lattice with even integer coefficients. 



The eight roots for SO(5) are given by 

E(50(5)) 



i^j = 1,2 



(12) 
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where is a unit vector in direction i. The short roots connect states along the coordinate axes while the 
long ones go diagonally. 

The same spinor also forms a representation of SO(4). However SO(4) has only the four long roots and 
one can see that the weights ±(i, i) do not talk to the pair ±(|, — i), which means the representation is 
reducible. We do not discuss it here. 

The simple roots are e! — e 2 and e 2 ; our root lattice is 2R = 2m(e 1 — e 2 ) + 2ne 2 = 2m'e 1 + 2n'e 2 , a 
square lattice of sides 2. A site on the decorated lattice is 2R + (±|, ±i). 

Since the spinor is self-conjugate, the original squares share corners with identical squares, and the resulting 



lattice is the square lattice. When links along the long roots are included as in Fig. 2(b) the structure is 
known variously as the square lattice with crossings (SLWC), checkerboard lattice, or planar pyrochlore. 
Note also that if you can hop from site j to site i on one unit cell, you can hop once more by the same 
amount to hit site j in the next unit cell. 



2.2.2. Sp(4) 

The weights of the 4-dimensional minuscule representation of Sp{A) are 

M = (±1,0),(0,±1). 



(13) 



Since the group is the collection of 4 x 4 symplectic matrices, this is also called the defining representation. 
(The same terminology applies, say to the 6-dimensional vector representation of 5*0(6), which is the group 
of 6 x 6 orthogonal matrices.) 
The roots of Sp(A) are 



E(Sp(4)) = ±e 4 ± e,- and±2e 4 i ^ j : 1,2 



(14) 



Now the long roots connect points parallel to the axes and short roots in diagonal directions, as shown in 



Fig. 3(a) Note that this just the rotated and rescaled version of 5*0(5) . 




Fig. 3. Planar pyrochlore and SP(A): The unit cell (a) for the defining representation of Sp(A) is the square rotated by 45 
degrees. The entire lattice (b) is the rotated version of the SO (5) spinor. 

This concludes the enumeration of lattices in d = 2. 
2.3. Lattices in d = 3 

With the warm up from d = 2 we can proceed rapidly through d = 3 where the candidates are 
SU(4)=SO(6), Sp(6), and 50(7). 



2.3.1. Sp(6) 

The only minuscule representation for Sp(6) is the defining 6-dimensional one. The weights are 

/x- (±1,0,0), (0, ±1,0), (0,0, ±1) 



(15) 



which form an octahedron which is self-conjugate. The decorated lattice is most readily visualized as corner 
sharing octahedra, which we will refer to as the octachlore lattice, depicted in Fig. [H The roots are 

E(5p(6)) = ±e, ± e 3 and ± 2e 4 i ^ j : 1, 2 or 3 (16) 

The short roots allow you to hop along the edges of the octahedron while the long roots take you straight 
across the unit cell to the antipodal point. This is just the d = 3 version of 5p(4). 




Fig. 4. The octachlore lattice of 5P(6): blue octagons are copies of the denning representation of Sp(6); red octagons are in 
the conjugate representation. Note that we have not drawn the bonds along the long roots. 

2.3.2. SO(7) 

Next we turn to SO(7). Only the spinor is minuscule. It has self-conjugate weights 

"=(444 (17) 

which lie at the corners of a unit cube. The roots 

Y,(SO(7)) = ±e; ± ej and ± i^: 1,2 or 3 (18) 

allow hops along edges, and diagonally across faces, but not along the body-diagonal. The root lattice L^r 
is cubic with edges of size 2. The decorated lattice is made of corner sharing unit cubes, with face diagonal 
hopping on every second cube as one proceeds along any of the three cubic axes (Fig. [5]) . This is one possible 
3 dimensional variant of the checkerboard lattice discussed in Sections 12.2.21 and I2.2.1[ we shall refer to 
it as the 3 dimensional checkerboard lattice. The other, which we will discuss in the next section, is the 
pyrochlorc. 

In higher dimensions the SO(2N + 1) spinor leads to the N-dimensional checkerboard lattice: an N 
dimensional cubic lattice with links on all face diagonals. 



2.3.3. SO(6) 

We obtain SO (6) from SO (7) if we drop the short roots ±e,: 

S(50(6)) = ±e 4 ± ej j : 1,2 or 3. (19) 

Consider the spinor representation of SO(7). Without the short roots ±ej, we can only flip the signs of 
the components of each weight two at a time. This means the 8-dimensional spinor of SO (7) breaks down 
into two irreducible representations with four weights each. The first has an even number of negative weights 

1 11 111 1 1 1 111 

' Xl= 2'2'2 '"* 2= 2'~2'2 " 3= 2'^ 2'~2 " 4= 2'2'~2 
and the other has these weights reversed and hence an odd number of negative weights. This is general: the 
irreducible spinor of SO(2N + 1) becomes two irreducible representations of SO(2N), called left and right 
handed spinors. 
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Fig. 5. The 3-d checkerboard SO(7) spinor lattice: blue and red cubes show the spinor representation and its conjugate (identical 
in this case). 

If you join the 4 points in either spinor multiplet, you will see the tetrahedra that form the weights of the 
SU(4) quark and antiquark representations. In other words the right and left handed spinors of 50(6) are 
the quark and anti-quark of SU(4). (This is why we will not study SU(4) separately.) The tetrahedra form 
the familiar pyrochlore lattice shown in Fig. [6l 




Fig. 6. The pyrochlore SO(6) spinor lattice: Blue tetrahedra are in the fundamental representation of SU (4), or the right-handed 
spinor representation of SO(6); red tetrahedra are in the conjugate representation. 

The extensions of this construction to higher dimensions yield corner sharing simplices that form the 
natural generalizations of the kagome and pyrochlore lattices [T5] . 

Since for SU(4) the difference of any two weights of the quark is a root and there are no others, the roots 
of SO(6) may just as well be written as 



a result we will invoke later. 

The third minuscule representation of S0(6) is the 6-dimensional defining representation with weights 



E(S0(6)) = ji< - Mj 1, 2, 3, or4] 



(21) 



/x = (±1,0,0), (0,±1,0), (0,0, ±1) 



(22) 
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just as in Sp(6). The decorated lattice again is made of corner-sharing octahedra, as shown in Fig. 0J 
However, without the long roots ±2ej, one can hop only along the edges but cannot jump directly from a 
point to its antipodal point. 

This concludes the enumeration of lattices. 

2.3.4. Structure factor 

Before we turn on the flux let us note that the structure factors for these lattices have a simple group 
theoretic interpretation. When we do a sum over all sites (indexed by I) we find 

5(k) = J2 eik r ' = X>' k ' 2R =^e ik m Tre ik H 

ieL 2R +fj, l 2R jefi l 2R (23) 

= SL 2R (k)x(k) 

where H; = Hi, the i th element of the Cartan subalgebra, whose j th eigenvalue is /z*-, and x(k) = Tre lkH 
is just the character of the representation. 

3. Flux Hamiltonians in d = 2. 

We will turn on fluxes by attaching arrows to each bond. The sense of the arrow will remain fixed as we 
move along bonds in any one direction. The hopping amplitude will be ±i if we go along (against) the arrow. 
Time reversal reverses every arrow, sending H — > — H . Since time reversal also reverses k, this means that 
in cases we study, 

ff(-k) = -H(k). (24) 

If all hopping amplitudes are pure imaginary, the flux through each triangle is ±7r/2. It is important to 
remember that the arrows themselves do not stand for physical quantities. For example a tetrahedron with 
uniform flux tt/2 coming out of each face is invariant under all symmetries of the tetrahedron although the 
arrows will look different if we say rotate the figure. This will, of course, be built into a PSG (Projective 
Symmetry Group) analysis [5]. 

3.1. SU(3) and kagome 

Let us begin with the first case, the kagome lattice, now with the flux shown in Fig[7] Note that as we 
go counter-clockwise around the (quark) triangles 1-2-3-1, we get a product (i) (+*)(— i) = +i. This is so for 
every quark triangle. The antiquark triangles will have the opposite flux. 




Fig. 7. Flux assignment on the kagome lattice. There is a phase factor of ±i as we move along (against) the arrow. 



Consider hops from the site numbered 1 in the central unit cell shown, to the sites numbered 2 to its left 
(on the same cell) and right (on the cell to the right). Hopping along (against) the arrows brings a factor 
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of ±i. We note that this state was one of several mean- field solutions found by [19] for uniform hopping on 
the kagome lattice. 

In terms of creation and destruction operators c and we have, in obvious notation, a contribution to H 



H 



1,21 



(25) 



2 (k) = 4a (ie*-^-^ + h.cj - (1 2) 
= (c|ci + c|c 2 ) sin(k ■ ai 2 ) 
where we have dropped an overall factor of 2, suppressed the k dependence of the operators, and as before, 

£*12 = Mi - JU 2 ( 26 ) 

Upon noting that c}c2 = Ea 12 , the generator corresponding to the root a.\2, we see that when all hops are 
included we get what we shall the canonical form 



H(k) = ^sin(k^ 



E f 



E. 



a 



(27) 



where the sum is over E + , the positive roots, which have the form — with i < j. 

In this paper we will often modify the canonical form in two ways: replace sin(k • a) by k • a and attach 
various signs sq = ±1 in front of each term so H assumes the more general form 



H(k, Sql) = Sa sin ( k ' a )( E a + E a)- 



(28) 



The role of these signs is to modify the relative phase of the hopping amplitudes on various bonds, and 
therefore the fluxes. 

Putting in explicit values, we get for the canonical case 

,1 





sin a; 



SIM 





sin(-(a; + v / 3y)) - sin(^(a; - V3y)) 



sin(-(x + V3y)) 
-sin(-(x- V3y)) 




(29) 



where x and y stand for k x and k y respectively. 
The determinant of this matrix 



\H\ 



-2 sin x sin 



i (x-V3yj) sin Q (x + VSyj^j 



(30) 



shows zeros along the lines x = and x — ±^/3y, which are precisely the directions of the weights (or their 
negatives)! The determinant also has a simple form in terms of the three positive roots: 



\H\ = —2 sin(k • a^) sin(k • 0:13) sin(k • a.23) (31) 

We do not know how this generalizes for SU(N). But we do know how to understand the lines of zeros 
as follows. 

The Hamiltonian has the form 



5^sin(k.( Mi -Mj))(Sa w +^ M ) 



(32) 



i<j 



Suppose we set k = fi 1 . (All points on the simplex are the same and we pick one that is easier to analyze.) 
The simplex has the property that 



a i = j 
b iy^j 



(33) 
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It follows that the argument of any sine in which fi 1 does not appear will vanish and the ones where it 
does will have the same value a — b. The resulting matrix, proportional to sin(a — b), has non-zero entries 
only in the first row or column. This leaves a 2 x 2 submatrix of zeros which kills the determinant: 



sin(a — b) 



3=2 



(E, 



OLu 



= 



(34) 



This will happen for all SU (N) because each raising or lowering operator in the fundamental representation 
has only one non-zero entry. (Geometrically this is equivalent to the statement that the N-simplex has no 
parallel edges). The components of the Hamiltonian which do not vanish for k = fi 1 are 



sin(a 



N 

&)X)(2Su 

3=2 



4) 



(35) 



But all of these are hoppings to or from site 1, so that all non-zero entries lie in either the first row or the 
first column of the matrix. The resultant (N — 1) x (N — 1) null submatrix ensures that the determinant is 
zero at this particular value of k. 

Let us now replace sinx by x, since none of the key features are lost and the algebra is more manageable, 
especially in the problem of diagonalization. So we will set 



H(k) 



\ 



The determinant of this matrix is 



\H\ 



1 



x(x - V3y)(x + \/3y) 



(36) 



(37) 



which shows zeros along the same lines x = and x = ±V3y as before. 

If we change the sign of the term multiplying any of the generators, we get the same lines of zeros. This 
is expected since we simply reverse the flux penetrating every triangle, which inverts the spectrum. If we 
reverse the signs of any two of them, it makes no difference even to the flux. The reader can check the lines 
of zeros are not altered by any change of sign. 

As explained earlier, this is a problem where H(k) = — H(— k), and the pair of points at ±k together 
produce E — > —E symmetry With the Fermi energy at zero we are at half-filling, the relevant filling for 
mean-field solutions of the Heisenberg model. 

The energies themselves are fairly complicated and not displayed here. It turns out two of them never 
vanish away from the origin and one of them produces all the zeros: if we move around the unit circle, it 
vanishes six times. Linearizing near these zeros will produce one-dimensional Dirac fermions that will control 
the low-energy physics. ( Near the origin all six Dirac excitations will get mixed up.) If all this were part 
of a mean-field calculation, we would be looking at this Dirac field minimally coupled to a gauge field if we 
wanted to consider fluctuations. One could ask if the mean-field solution remains stable in their presence. 
These questions will be considered separately [2"U] . 



3.2. Flux Hamiltonians for 5p(4)=50(5) 



Recall that the groups Sp(4) and 50(5) are the same. The roots of one are the rotated and rescaled 
versions of the other and no new physics will come from looking at both. We will only work with 50(5) 
since it may be more familiar to the reader. 

Before writing down the hopping matrix we need to define the basis. The states are numbered 1 through 
4 in Fig. [8] with 1 = (|, 3) etc. We use a tensor product of two Pauli matrices, a and r, to operate on the 
two labels. We take as generators 
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^ei±e 2 — T [^61)^62] = = ' = 2 Cr " 1 



(38) 



3.2.1. The canonical Hamiltonian 

In this basis we choose the canonical form 



H(k) = ^(k- a) (£«+#, 






y 


X 


se + y 


y 





y-x 


X 


X 


y-x 





-2/ 




- 2/ z 


-2/ 






(39) 



/ 

The orientation of the arrows corresponding to this H are shown in Fig. [8] In other words, rather than write 
down some arrows and deduce the hopping matrix from these, we are writing down a canonical matrix in 
the Lie algebra and asking what hopping elements it implies. 




Fig. 8. Flux assignment for the canonical SO (5) Hamiltonian. There is a phase factor of ±i as we move along (against) the 
arrow. The flux in each triangle alternates as we go counterclockwise in the case depicted. 

If we look at the flux in each triangle we find that it alternates: triangles sharing a face diagonal, such 
as triangles 132 and 342 in Fig. [8l have the same flux, while triangles not sharing a face diagonal (triangles 
132 and 142) have opposite flux. The determinant here is 

\H\=4(x 4 -x 2 y 2 + y 4 ) (40) 

It has permutation symmetry but not full rotational symmetry. It has no zeros anywhere away from the 
origin. We do not know a simple way based on group theory to understand this. 
The characteristic polynomial is 

£ 4 -4£ 2 (x 2 + y 2 )+4(x 4 -ry + y 4 ) (41) 

resulting in the particle-hole symmetric spectrum: 

E = ±V2\Jx 2 +y 2 ± V3xy (42) 
The reason behind the symmetry E — > — E is the matrix 
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G 






-1 




1 

-10 



A 











(43) 



with G 2 



I. Since 



G ■ H ■ G = —H 



(44) 

it follows that H and —H have the same spectrum. 

The anatomy of the operator G is interesting. Suppose we wanted to manufacture an operator that reversed 
the sign of H by conjugation. We could accomplish this by a parity operation that exchanges each weight 
with its negative- this should flip every hopping term. However there are many ways to flip the weights since 
we can take each state to its parity reversed state, times any unimodular phase factor, which must be a sign 
if we want the hoppings amplitude to be ±«. Suppose we picked a G' with all positive signs: 



G' = 







10 
10 
10 



(45) 



We would find 



G'HG'- 1 = 



(° 


-y 


X 


X- 


ry 


-y 





y — x 


X 




X 


y-x 





V 




\x + y 


X 


y 





) 



(46) 



This is clearly not —H. However, if we go back to the lattice and ask what the corresponding hopping 
amplitudes are we will find that all the fluxes are reversed, though some arrows are reversed and some are 
not. Since —H also has all fluxes reversed (reversing every bond will reverse the product over every triangle) 
the two must be gauge equivalent. It turns out that appending minus signs to states 3 and 4 is one way 
to flip the arrows that needed to be flipped. The operator G is the product of G' and a diagonal matrix 
that multiplies 3 and 4 by minus signs. This is to be expected in a gauge theory, where the symmetry is 
projectively realized [5]. 

Since such a procedure will work for any self-conjugate representation, we will not explicitly construct the 
operator in future occasions. 

We can understand now why the spectrum of H had the full set of lattice symmetries even though the 
flux alternated. Under any of the lattice symmetry operations, we either left the flux alone or reversed it. 
Neither affects the determinant since these operations at worst exchange E — > —E, which has no effect on 
the spectrum. This feature will be seen again when we consider other groups. 

Since H is 4 x 4, and the characteristic equation is even in H , it satisfies an equation of the form 



(H 2 -f(k)-lY=g(k)-I 
where / and g are scalar functions. Eq. (|47| reduces to the Dirac form if g(k) — and / is constant. 



(47) 



3.2.2. The non-canonical Hamiltonian 

Let us now change the sign in front of any of the terms in Eq. 1391 It turns out that the determinant is 
sensitive only the relative sign of the two long roots that reach diagonally across the square. Here is what 
we get when we flip the coefficient of the E corresponding to the root ei — e-2 : 
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o 
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x 



x + y 



V -(y - x) x 

x -(y-x) —y 



l x + y x 



-I) 







(48) 



If we compute the flux now, we find it is uniform in all triangles, hence the subscript in H un i . This in turn 
means that H un i will be invariant (up to gauge transformations) under symmetry operations of the lattice. 
The characteristic polynomial is 



E 4 



AE 2 x 2 



4E 2 y 2 



4x 2 y 2 



(49) 



resulting in the spectrum 



E = ±\/2 V x 2 +y 2 ± a/z 4 + x 2 y 2 + y 4 (50) 
which has E — > —E symmetry because G once again anticommutes with H. The spectrum has zeros along 
the lines x = and y = 0. Note that along these directions, ay • e x = 1 or 0. As in the kagome case, the lines 
of zeros are axes of symmetry of the unit cell. All edges not orthogonal to these axes have equal projections 
onto them (up to sign). With signs as in H un i, this results in a pair of identical rows in the Hamiltonian - 
and thus lines of zeros along the short roots for both the linearized and lattice versions of H. 



3.2.3. The Hamiltonian with unequal hopping 

In the above we have considered a more general case 



H(k, 8a) 



s a sin(k • a)(E a 



El) 



(51) 



where sot = ±1 is a possible sign. In this problem there were essentially just two choices, the ones with 
uniform and alternating fluxes, determined by the relative sign of the two long hops. 

In a problem like 5*0(5) = Sp(i), where there are roots (i.e. bonds) of two different lengths, we could also 
play with the relative strengths of the hopping across long and short bonds. There is no obvious inspiration 
from group theory on how to choose from the continuum of possibilities, though only some choices will give 
SU(N) mean-field solutions. The only consolation is that only two different lengths are allowed for the roots 
of any semi-simple Lie algebra and among the cases we study this happens only for SO(2N + 1) and Sp(2N). 

We just mention one extreme case where long hops are set equal to zero: 



H 



short 



(k) 



f° 


.'/ 


X 





y 








X 


X 








-y 


1° 


X 


-y 






(52) 



Now the determinant is 



\H short \ = (x 2 +y 2 ) 2 (53) 

which describes two Dirac points at the origin. Indeed, this is just the flux phase on the square lattice 
originally described by 0]. The unit cell is of course twice as big as it needs to be, so that the two Dirac 
points of the traditional unit cell have both come to the origin. 

More generally, with a magnitude c for the long hops and relative signs all positive, the eigenvalues are 

E = ±\J{1 + c 2 ){x 2 + y 2 ) ± 2cxy\/c T +^ (54) 
and the spectrum has one Fermi point at the origin. For relative signs chosen as in H un i , the eigenvalues are 



E = ±\J (1 + c 2 )(.t 2 + y 2 ) ± 2c v / x i + y i + c?x 2 y 2 



(55) 



1G 



At the special values c = 0, y/2, these give a Dirac spectrum. For all c ^ 1, the Fermi surface is a single 
point at the origin. Values of c corresponding to mean-field solutions are given in Section [5.11 



4. Flux Hamiltonians in d — 3 



Luckily we have to consider just three groups: 50(6)=5J7(4), 5P(6), and 50(7). The only minuscule 
representation of the latter is the spinor. There are three minuscule representations for 50(6): two spinors 
(quark and antiquark of 5C/(4)) and the six-dimensional vector representation. 5P(6) has one minuscule 
representation- the defining one. 



4.1. Flux on the 50(6) spinor lattice. 



The weights forming the tetrahedron are 

111 111 111 111 

^ 2'2'2 ^ 2 _ 2~' _ 2'2 2' _ 2' _ 2 ^ 4 2'2' _ 2 
The positive roots are, in terms of orthogonal unit vectors, 

ei±ej j>i = l, 2, 3. (57) 

However, as pointed out in Eq. (|2ip. it is more convenient to note that since this is also an SU(4) quark 
representation, we could write them in terms of the weights (|56[) as 



In view of what we saw in d — 2 we are going to admit the more general case 

#(k, s a ) = ^2s a sin(k • a)(E a + E j a ) 

s+ 

where sa = ±1 is a possible sign in front of each term. 



(58) 
(59) 



4.1.1. The canonical Hamiltonian 

If we pick all signs positive, (which means the arrow always goes from a site with a lower index to one 



with a higher index) we obtain the flux assignment in Fig. 9(a) This gives the canonical Hamiltonian 



H 








x + 
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+ z 
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+ z 


x - 


Yy 







z 


— X 
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- y 


1H 


- z 
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X 
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Yz 
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with determinant 



\H\ = Ax 4 - 4xV 



4x 2 z 2 + 4y 2 z 2 



(60) 



(61) 



Note that it lacks the discrete symmetries of the lattice. This is to be expected since the flux on each face 
is not the same. 

Consider its zeros. It vanishes along the weight directions, k oc [i^ This is to be expected since the right- 
handed spinor is also the SU(4) quark representation and we have seen that for SU(N), because the weights 
form a simplex, when k cx only terms corresponding to roots involving ii i remain (and that all have the 
same coefficient in front). The rest vanish, so that H has just one nonzero row or column. 

But we find in addition that there are entire planes along which there are zeros. For example for any 
linear combination k = afi 1 + bfi 2 or k = a\x x + 6^t 4 , the determinant vanishes. However it does not vanish 
for k = afi 1 + b/ji 3 unless a = or b = 0. This variation is to be expected since the flux is not symmetric on 
the tetrahedron. 
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Once again if we can use more powerful group theoretic methods to know when determinants of certain 
elements of the Lie algebra of the type of Eq. (|59|) will vanish, we will be able to anticipate this result rather 
than just observe it. 





(a) 



(b) 



Fig. 9. Possible flux assignments to the pyrochlore. (a) Flux assignment breaking rotational symmetries, (b) Flux assignment 
preserving rotational symmetries. 



4.1.2. The uniform monopole case 

If we flip the sign of the 24 and 42 matrix elements (corresponding to the root /x 2 
the flux assignment of Fig. |9(b)| The Hamiltonian is 



/x 4 = CK24) we obtain 



x + y 
x + y 
y + z z — x 



y + z x + z 
z-x -(z-y) 
x — y 



1 x + z -(z -y) x - y 







(62) 



with determinant 



\H mono \ = Ax A + 4y 4 + 4z 4 - Ax 2 y 2 - Ax 2 z 2 - Ay 2 z 2 



(63) 



which has the discrete symmetries of the lattice. The subscript on H mono reflects the fact that the flux is 
the same on all faces of the tetrahedron and comes from a unit monopole at its center. 

The reader may well ask how many more such signs are we going to play with. Luckily we are done. 

To understand this, we need to transcribe the Hamiltonian to the corresponding factors of ±i on the edges 
of the tetrahedron. As mentioned above, the case with all sex = 1 has a factor of i if we move from a corner 
to another with larger index (and a — i if we move the other way). It is readily verified that the two faces 
not involving the bond 24 have an outward flux of ^ir (a factor i around the triangular faces) and the other 
two the reverse. Clearly the choice of matrix elements violates the discrete symmetries of the tetrahedron. 

On the other hand if we flip the coefficient of the 0:24 term, we get an arrangement with all outward fluxes 
equal to it/2 and we are led to H monoi with a tetrahedrally symmetric determinant. 

Other choices of sign will only yield one of two options: the flux is uniform (could be ±^7r) and the 
determinant is symmetric, or the flux assignment breaks the symmetry with two positive and two negative 
faces. Different choices for the latter will correspond to determinants in which the asymmetric roles of x, y 
and z are interchanged. 

We will now elaborate further on the case Hmono, which describes uniform flux, as it has various nice 
properties. It has been discussed elsewhere as an interesting mean-field Hamiltonian for the SU{2) Heisenberg 
model on the pyrochlore lattice [21122) . To make contact with existing literature on this problem, we will 
briefly revert to the custom of referring to momentum components as k x or k v rather than simply x or y. 

The energy levels of H mono are 
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E{k) = ± 



(64) 



(i<j) 



The spectrum has E — > —E symmetry since there is a matrix G obeying G 2 = I that anticommutes with 

HrnonoG — GHjnono (^^) 



Emono • 



where 



G = 



v/3 



111 



1 
-1 



-1 1 -1 







(66) 



That a matrix G which ensures E — > — E symmetry should occur is less obvious than in the 50(5) case, 
since tetrahedron is not inversion symmetric (self-conjugate). The inversion operation maps the tetrahedron 
formed by the right-handed spinor representation of 50(6) to that formed by the left-handed representation. 
Hence G is not a simple geometric operation on sites in the unit cell. 



We are tempted to cast H mono in Dirac form 

Emono — ®-xkx 



atyky + a z k z 



(67) 



ijk 



since G seems to be like the matrix (3 which anticommutes with the three a's in the Dirac Hamiltonian. 
However the resemblance to the Dirac case is not complete because a's do not form a Clifford algebra and 
Hmono i s n °t a multiple of the unit matrix. 
What one finds is 

[oti, aj] + = 25ij + \/3\t 

In other words, the anticommutator of the a's is proportional to the unit matrix plus some amount of W's, 
and the W's obey a Pauli algebra. Thus if we square H(a), move the stuff proportional to the unit matrix 
to the left hand side and square again, we will end up with a multiple of the unit matrix. Indeed this is so: 



(68) 



(H 2 - 2k 2 ) 2 = 12(k 2 x k 2 y + k 2 x k\ + k 2 y k 2 z ) 

where the subscript on H and the identity / have been suppressed. 

That we should end up with the form encountered in the d = 2 50(5) case of Eq. (|4"7| 

(H 2 f(k) ■ I) 2 = g(k)I 



(69) 



(70) 



is due to the same reasons: the characteristic polynomial P{H) is even and of fourth order in H, i.e. quadratic 
in H 2 . It can therefore be cast in the form Eq. ([70)) . To get all details of / and g we would of course need 
to actually evaluate P{H): 

H A - 4H 2 (kl + k 2 + k 2 z ) + 4(fc4 + ky + kt- k 2 x k 2 - k 2 k 2 z - k 2 z k 2 x )) 



P(H) 



0. 



(71) 



The anticommutator algebra in Eqs. (|6"8]) stems from the fact that the Hamiltonian lives in the Lie algebra 
of the right handed spinors of 50(6), which is also the quark of SU(4). 

Recall that it is possible to write the generators of SO(N) in the spinor case in terms of the Dirac 7- 
matrices: cr^v, which generates rotations in the \x — v plane, may be expressed as — ^7^7^. Although the 
7 matrices are 8x8, bilinears in them like a^ u form reducible representations with two 4x4 blocks, these 
being the quark and antiquark of SU(4). The two blocks are eigenstates of 77 = 171 • ■ - 76 with eigenvalue 
±1. If we want the quark we can work with these 8x8 matrices and focus on just the top left hand corner. 
In this block 77 is just a number equal to 1. 

Consider the following operator 



H = 171(76 - n)k x + 173(72 + le)ky + 175(74 + 72)^2 = Hr®H l 



(72) 
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where Hr and Hl are 4x4 blocks corresponding to right and left handed spinors, or quark and antiquark 
representations . 

With a judicious choice of basis for the 7 matrices its upper left-hand corner, Hr is just our Hamiltonian 
ot x k x + OLyky + a z k z . Thus if we do not stray from this block we can view the a's as bilinears of 7 matrices. 
Not so obvious is the fact that the Ws which come from two powers of a are also bilinears in 7. 

The closure under anticommutation of the <j^ v or the a's and Ws is a special property of SO (6). In 
general, if you multiply two of them you will get something quartic in the 7's even after some of them reduce 
to quadratic terms upon invoking = I. The quartic ones can be rewritten as 77 times a quadratic, upon 
inserting the square of the "missing" two 7 matrices. In the sector with 77 = 1, these are just quadratic in 
the 7's. 

4.2. Flux on the SO(6) vector lattice. 

In the defining vector representation the generators are represented as follows in terms of canonical creation 
and destruction operators c and c': 



H, 



= c\ci — „ c_ 



E. 



i = 1,2,3 



c\.-c-i i < j < 3 



(73) 



with generators of negative roots defined as the adjoints of the positive ones above. 

4.2.1. The canonical Hamiltonian 

In this basis the usual sum over positive roots with all coefficients positive yields the matrix 

A 







H = 






x — y —x — y 
x + y y — x 
x — z —x — z y — z 
\^x + z z — x y + z 



x—y x+y x—z x+z 
—x — y y — x —x — z z — x 
y — z y + z 

— y — z z — y 

-y-zO 
z-y 



(74) 



where the rows and columns are numbered as follows: (1, 0, 0, ) = 1, (—1, 0, 0) = —1, (0, 1, 0) = 2, ...(0, 0,-1) = 
—3, the components being just the eigenvalues of Hi, H2 and H3. The site labels and corresponding factors 
of ±i are shown in Fig. [101 Note that the flux alternates from one face to the next. 




Fig. 10. Flux assignment to the octachlore in accordance with signs of the group generators of the vector representation of 
,50(6). 
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The determinant vanishes identically because there are two zeros at every k. If we pull them out we find 

\H\ = 48{x 2 y 2 + y 2 z 2 + z 2 x 2 ) *0*0 (75) 

that is to say, the product of the nonzero energies is 48(x 2 y 2 + y 2 z 2 + z 2 x 2 ). 

Why does H have all the discrete symmetries when the flux alternates? The answer is that any rotation 
is equivalent to a change of the sign of the overall flux, which in turn corresponds to time-reversal, and does 
not affect the determinant in a problem with E — > — E symmetry. 

Extra zero-energy bands occur when any two coordinates vanish, i.e., along the axes, which corresponds 
to the direction of the weights. We can understand this to the extent we could understand the SU(N) and 
SO (5) cases. If 

H = J> • (e, ± Bj)) [E e%±ej + El i±e .] (76) 

it follows that if we set k = ei say, only roots of the form ei ± will survive and that too with the same 
coefficient. The matrix will have only two non-zero rows and columns - for the sites at ±ei in the unit cell. 
With the flux assignment of (|74[). one row is exactly the negative of the other, resulting in two extra zero 
energy bands in both the linearized and the lattice Hamiltonian. 

Near any line of zeros we can define a 2 dimensional Dirac field, except near the origin when they all 
collide and modify each other. 

Again there is a matrix G which anticommutes with the Hamiltonian, and acts upon the unit cell as the 
inversion. Its existence results from the fact that the unit cell is inversion symmetric, while the directions 
of all fluxes are reversed by inversion. 

4.2.2. The non-canonical Hamiltonians 

We could append signs for each term, but this gives spectra which break the lattice symmetries. We have 
not looked deeply into what kind of zeros result in that case. 

We did however note the following. Suppose we start with a hopping problem on an octahedron with 
uniform flux in every face as in Fig. 1111 




Fig. 11. Flux assignment to the octachlore preserving rotational symmetries. 

When we extracted the H(k) for that problem we found it could not be written in terms of generators 
of SO(6). It is important to understand why we have this problem here but did not when we considered 
the spinor of 50(6) = quark of SU(4). There each root or generator connected only two states. If we did 
not like the sign of the matrix element given by group theory we just put a negative sign in front of that 
generator using sa. But here, each root connects two pairs of points, corresponding to parallel edges of the 
octahedron. For example 

E ei+e2 = c\c- 2 - 4 c -i ( 77 ) 
connects points labeled (—2, 1) and (—1,2) in Fig. [TT] with opposite matrix elements. If we do not like the 
relative sign, we cannot do anything about it. This is exactly what happens in the case of the octahedron 
with uniform flux in every face. To describe it, we would have to use the generators of the much larger group 
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SU(6). But that is not the game we are playing: we want to work within a group, SO (6) being the operative 
one here. 

This could have happened to the 50(5) spinor, whose short roots connected opposite sides of the square 
with same sign for the horizontal roots ±ei and opposite signs for the vertical roots ±e 2 . Luckily this choice 
of signs corresponded to the case of interest. 



4.3. Sp(6) 

The defining representation of Sp(6) is the same octahedron as in SO (6) with the same weights. The 
generators can be written in terms of creation and destruction operators as 



E, 



c\ci 
cjc 



cL.-c_ 



E 2e% = 2c] c 



i = 1,2,3 

: i < i < 3 



(78) 



with negative roots being given by adjoints of the above. 



4.3.1. The non-canonical Hamiltonian 

It is interesting to consider first the canonical H with all signs positive and only the short roots. (As noted 
before when we have two different root lengths, we have the freedom to chose the scale of each type of term. 
Keeping only short roots is an extreme case.) We find 



H 



short 



(k) 



/ x — yx + yx — zx + z 

x + yy — xx + zz — x 

x — y x + y 

x + yy — x 

x — zx + zy — zy + zQ 

yX + z z — x y + z z — y 



y - z y + z 
y+z z-y 





(79) 



The determinant has the value 



\H short \ = -32(x 2 + y 2 )(x 2 + z 2 )(y 2 + z 2 ) 



(80) 



that is to say, zeros along the weights. The logic is the same as in SO(6) since the long roots that distinguish 
between them have been suppressed. Note however that matrix elements are different now: the two pairs of 
states connected by a generator do not always have opposite matrix elements. Thus in this case along the 
axes there are 2, rather than 4, zero energy bands. 

What is surprising is that the energies do not change for any choice of signs ! 



4.3.2. A non-canonical Hamiltonian with unequal coefficients 
Consider the following matrix involving the long roots: 
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(81) 



y —x — z z — x y + z z — y 2z J 

where the subscripts remind us of two ways in which it differs from the canonical form: the x + z term has 
a minus sign relative to the canonical form, and the hopping matrix clement for the long roots is half as big 
as the canonical one. As for the latter point, consider the term 2x. It indeed equals k • 2ei, but the generator 
E% ei = 2c{c_i has another two in it. So this term should have been Ax. But with the choice of sign and 
hopping (|5Tj) we get 

\H W ,±\ = -l&{x + yf(x + zf{y + z) 2 (82) 
which has zeros along planes x + y = etc. 

When we put in the canonical strength [Ax etc.) we did not find any interesting spectra for many choices 
of sign that we tried. 



4.4. SO(7) spinor 



Recall that the only minuscule representation of SO(7) is the spinor and that the lattice we associate with 



it is cubic, with face diagonals but no body diagonals (Fig. 12(a) ). The generators in this representation can 



be expressed in the direct product space of three Pauli matrices: 

E± ei = ier 3 <g> r 3 <g> a± E± e2 = ier 3 ® T± ® 1 

I 1 

E± es = -<r± 1 1 E ei±e2 = ±-l (g) t + (g) a± (83) 

1 I 
E ei±e3 = ±2 a + ®T 3 <g>a± E e2±e3 = ±2 CT + ® t± ® 1 

where the labels 1, 2, and 3 correspond to the directions of the co- ordinate axes. 
The matrix 

G = (T 2 ® Tl ® «2 (84) 
acts as an inversion operator on the unit cell and anti-commutes with all of the symmetric generators 
E a + E- a so that the spectrum has symmetry under E — * — _B no matter what the signs. 



4.4.1. The canonical Hamiltonian 

If we ask what hopping amplitudes are associated with the canonical case we find that each square 
plaquette has ir flux. 

The generators Ei±j fix the flux through the triangular plaquettes to be ±7r/2. The form of these generators 
dictates that two opposing pairs of triangular faces will have diagonals with the same orientation; the third 
will have diagonals with opposite orientations. It turns out that in this case a uniform flux through the 
triangular plaquettes is impossible. 

If we choose all signs to be positive, then three of the cube's faces have flux ir/2 outwards through all 
triangular plaquettes, and the remaining three to have flux —tt/2. 

The zeros of energy can be found from 

l^l = ^rS x2 + y 2 + z2 - 2 ( xz -y z ~ ^)) 2 ( x -v + z ) A ( 85 ) 

25o 

Here not all cubic symmetries are preserved, but permutations of the x, —y, and z axes (corresponding 
to rotations of the cube about the (1,-1,1) body diagonal) map positive fluxes to positive fluxes and vice 
versa, so that some of the cubic symmetries are preserved. 
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4.4.2. The non-canonical Hamiltonian with alternating flux 

As in the octahedral case, the other interesting case is the alternating flux pattern shown in Fig. 12(b) 
in which rotations by 7r/4 about the x, y 7 and z axes reverse the signs of all fluxes. In this case, we append 
minus signs to the x + z, x — y, and y — z terms. The energies are remarkably simple: 

E = ±^V3{x±y±z) (86) 

which vanish along the planes x = ±y ± z. For example, if the sites of the unit cell are labeled 1 — 8 as 
shown in Fig. |12(b)| momenta in the plane x = — y — z obey k ■ ot.\j = — k • c*8j, producing a pair of linearly 
dependent rows in the Hamiltonian. In the alternating flux case the cubic symmetries are preserved since a 
7r/4 rotation reverses the signs of all fluxes, which is gauge equivalent to reversing the signs of all hoppings 
and thus the sign of H. Invariance under E — ► —E thus ensures that this is a symmetry. 

These two possibilities are the only ones preserving the permutation symmetry of the x, y and z axes. 



(a) (b) 

Fig. 12. Two possible flux assignments for on the SO(7) lattice. Each square plaquette has flux it, while the triangular plaquettes 
on each face may be chosen to have equal (a) or alternating (b) flux. 



5. Comments 



5.1. Relevance to mean-field solutions of the Heisenberg model 

We will now revisit the Hamiltonians of Sections [3] and 0] with a view to asking whether they are, in 
fact, mean-field solutions of the SU(N) Heisenberg model. A few of these have been discussed previously 
in the literature. The kagome Hamiltonian of Section l3Tl was identified as a mean-field solution by [19] . 
Several authors have discussed SU(N) mean-field states on the square lattice with second neighbour hopping 
9 8,23 24 , though these have focused on the gapped chiral spin state quite unlike that of Section l3~2{ we 
believe that previous mean-field studies on the checkerboard lattice [35] have been restricted to dimerized 
states. 

We wish to extend these results and establish that many of the Hamiltonians discussed above are mean- 
field solutions. To show this, we must argue that the mean- field equations admit solutions in which the 
hoppings t a are purely imaginary, and hoppings along roots of equal length are equal in magnitude. 

The mean-field equation can be written [5] 

£ p^ H * (k)e ~' k ' r " (87) 

JqelBZ B („) (k)<0 I V )\ 

where Jy is the spin-spin coupling between sites i and j for the original Heisenberg Hamiltonian, and is 
the hopping matrix element between these sites for the mean-field Hamiltonian. If the spectrum is invariant 
under k — > — k, then H(k) — —H(~k) guarantees that the real part of the integral vanishes. Further, if 
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the spectrum preserves the lattice symmetries, then the integral will clearly have the same magnitude for 
all roots onj of a given length. In cases with two different length roots, the mean-field equations specify a 
particular relative hopping strength. 

The reader should be warned, however, that some of the Hamiltonians we have studied have large manifolds 
of zero energy, such as planes of zeros. These may lead to a divergent integral in (|87p . in which case the 
Hamiltonian is not a mean-field solution. 

Such exceptions aside, many of the flux configurations discussed in Sections [3] and [5] do clearly give 
self-consistent mean-field Hamiltonians. Consider first the lattices in which all links have the same length 
(and are related by lattice symmetries): the kagome, pyrochlore, and octachlore. In these cases (|87| gives 
equal hopping amplitudes on all links provided that the spectrum does not break the lattice symmetries. 
The kagome Hamiltonian has been discussed by 19J; it is somewhat exceptional among the Hamiltonians 
we consider in that its spectrum is not inversion symmetric at fixed k. The monopole Hamiltonian of 
the pyrochlore lattice ( Section |4~L2| . as well as 50(6) (Section |4~2~T) and short-root 5P(6) (Section |4XT|) 
Hamiltonians on the octachlore lattice have spectra that preserve lattice symmetries and thus naively should 
be mean-field states by the argument above. We have verified numerically that the monopole and short- 
root 5P(6) Hamiltonians are viable mean-field solutions. The 5*0(6) Hamiltonian on the octachlore lattice, 
however, has two bands of zero energy, causing (|87[) to diverge; it is not a mean-field solution. 

The lattices constructed from the spinor representation of SO(2N + 1), as well as the canonical lattice 
of SP(2N), have roots of two different lengths. In all three examples discussed here, Hamiltonians with 
hopping only along the short roots are mean-field solutions; for 50(5) this gives the flux state of [3] and for 
50(7) the 3-d version thereof. If we are interested in Hamiltonians with non-zero hopping along the long 
roots, we must find the ratio ti/t s consistent with Eq. I|87[) . 

A summary of the allowed values of ti/t s is given in Table IBTT1 For 5P(6) generic values of ti/t s do not give 
symmetric spectra, and we find no mean-field Hamiltonians with t; > 0. The alternating flux Hamiltonian 
of 50(7) has a symmetric spectrum for general ti/t s ; however for many ti/t s two-dimensional surfaces of 
zero energy cause the integral to diverge and we are unable to find a mean-field solution with ti > 0. Both 
flux assignments discussed in Section 13.21 for the 50(5) spinor lattice give symmetry-preserving spectra 
whose mean- field U /t s can be calculated numerically. The ratio i; jt s at mean- field depends on the relative 
magnitudes of the spin-spin coupling J//J s in the original Heisenberg Hamiltonian; solutions with ti > 
exist only for sufficiently large Ji , as shown in Tabic 15.11 





SO(5) (alt) 


50(5) (uni) 


SO (7) (alt) 


5P(6) 


Jl/Js 


ti/t a 
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ti/t s 
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0, 1.59 


0,2.28 
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.371 


1.66 








.8 





1.28 








.7 





1.07 









Table 1 

Relative strengths of hopping along the long (<;) and short (t s ) roots as determined by the mean-field equations for the SO (5), 
SO(7), and SP(6) hopping problems. For SP(6) and SO(7) we find only the t; = solution at mean-field level. For SO(5) we 
find consistent mean-field solutions with i; > for Ji/J s > .682 in the uniform flux case, and .891 in the alternating flux case. 

Needless to say, this analysis does not preclude the existence of other flux assignments leading to a lower 
mean-field energy. Indeed states with lines and planes of zeros, as many of our examples have, are often 
energetically disfavoured at mean-field [19] due to the large phase space near E — relative to gapped 
or mostly gapped states. Also one must bear in mind that dimerized mean-field states of lower energy 
inevitably exist |ll|12j . However, as pointed out in Section [LT| corrections to the mean-field solution for 
N < oo often alter the relative stability of various mean-field states, so we should not take this issue too 
seriously. Among the mean-field solutions discussed here, an interesting example of this is the monopole 
Hamiltonian of the pyrochlore lattice. It corresponds to the lowest energy symmetry preserving mean-field 
solution to the SU(N) Heisenberg model, and has lower energy than the dimerized mean- field ground states 
after Gutzwiller projection is used to enforce the constraint of single occupancy |22j . 
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5.2. Hamiltonians beyond the linear approximation 

This work has focused on the linearized versions of lattice Hamiltonians, in which we have replaced 

sin(k • r) -> k • r. (88) 

However, many of the interesting properties of the spectra are unaffected by this substitution. 

First, surfaces of zero energy which are related to symmetries of the unit cell will not be affected. Recall 
that we find several zero-energy surfaces along directions of the unit cell for which k-r^ takes on values ±c, 
for some constant c. If k • is replaced by sin(k • ry), the only effect on the Hamiltonian at these points 
is to change the value of the constant c; hence the zero eigenvalues remain. As discussed above, this yields 
lines of nodes along the weight vectors for the SU(N) lattice model for any N. The zero-energy manifolds 
of the alternating flux 50(5) and 50(7) Hamiltonians, and the 50(6) and SP(6) Hamiltonians discussed 
above for the octachlore lattice, are also preserved under In other words, all zero-energy surfaces listed 
in Sections [3] and 0] which reflect symmetries of the lattice unit cell are unaffected by the substitution ([85]) . 

Second, the symmetry of the spectrum will remain. Symmetries in the spectra occur when flux is assigned 
in a way that preserves the lattice symmetries, and the substitution (|88p cannot alter the symmetry properties 
of the state. 

Finally, it is interesting to note that on the lattices with inversion-symmetric unit cells (namely lattices 
related to representations of SO(N) or SP(2N)), the operator G which anti-commutes with the linearized 
Hamiltonian also anti-commutes with the lattice Hamiltonian. This happens because G in these cases is 
simply the inversion operator multiplied by an appropriate gauge transformation, and inversion maps every 
edge to another edge associated with the same symmetric generator of the Lie group representation. Thus G 
anti-commutes separately with all of the symmetric generators - and hence also with the lattice Hamiltonian 
J2 s+ (E a + El) sm(k-r a ). 

5.3. Extensions to Higher Dimensions 

We have already noted that the generalization of our Hamiltonians to d > 4 is problematic. For com- 
pleteness we note here that the lattice construction described in Section [2] can be applied to the appropriate 
representations of the Lie groups discussed above in arbitrary dimension. Assigning a hopping of ±i to 
each directed edge will result in a Hamiltonian related to the group generators by Eq. (|28p . for which 
H(k) — — H{— k). Additionally, for all of the cubic lattices (SO(N) spinor, SO(2N) vector, and the defining 
representation of Sp{2N) ) a matrix G can be found which anti-commutes with the Hamiltonian, leading to 
a time-reversal invariant spectrum. 

In general, however, the symmetry operations of the resulting unit cell make it impossible to assign flux in 
such a way that the lattice symmetries are unbroken. The notable exception is the SO(2N + 1) spinor case 
with only the short roots, where all fluxes are tt = — 7r. This gives the N dimensional Dirac Hamiltonian. 

6. Conclusions and Outlook 

In this paper we constructed a class of lattices inspired by the root and weight systems of Lie algebras. 
The lattices had as their unit cells minuscule representations of the standard Lie groups which decorated the 
underlying lattice L>2R, elements of the root lattice with even coefficients. We observed that the lattices could 
equally well be viewed as decorations of L2R by the conjugate representation, which shared corners with the 
original one. While construction works for any rank r we stuck to r — 2, 3 since these were experimentally 
accessible and because these allowed an unambiguous assignment of flux on the triangular faces of the unit 
cells. Remarkably, they also correspond in many cases to known lattices like the pyrochlore, kagome or 
checkerboard. Even our octachlore lattice is a motif in the perovskite structure. 

We find this last aspect enticing, for it hints that it may be possible to relate more physics on these 
lattices to the underlying Lie algebras. Indeed, as we were finishing up this work we came across recent 
work by Arovas [26] who constructs generalized AKLT models on the kagome and pyrochlore lattices which 
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naturally involve local degrees of freedom that live in the fundamental representations of 5(7(3) and SU(4) 
respectively. 

However, our own work makes a different connection. We considered hopping Hamiltonians which, when 
written in momentum space, were elements of the Lie algebra, linear in the momentum k for small k, and 
obeyed HQs) = -#(-k)0D By varying the signs in front of each generator we could alter the fluxes in the 
faces of the unit cell. We found Dirac or Dirac-like spectra at points, lines and even sheets. The locus of 
the zeros had strong ties to the directions of the weights or roots. We could anticipate and thus understand 
some of them using ideas from Lie algebras but often were just able to draw attention to them. It seems very 
likely that an assault using ideas from Lie algebras can yield further understanding. To begin with one must 
employ a more systematic way to represent weights and roots in dual bases: simple weights for the former 
and simple roots for the latter. One should also use color groups to classify symmetries of this problem 
where the triangular faces of the unit cell are colored with flux ±7r/2. Of all the properties associated with 
H, the determinant seems most likely to yield to group theoretic methods. It has uniformly proven to be a 
much simpler and more symmetric function of the momenta than individual eigenvalues. 

The spectra often had E — > —E symmetry. For self-conjugate representations we could fully understand 
this feature and indeed use our understanding to construct an operator G that anticommuted with the 
Hamiltonian and explained this feature. 

While such hopping problems typically arise as lattice regulators for continuum theories or as mean-field 
theories for quantum spin models, in this paper we have studied them in their own right. While we did 
observe that most of them are candidates for interesting mean-field theories of quantum Heisenberg models 
on the same lattices, a fuller investigation of the fluctuations would be required to establish their value in 
that setting. 

We hope that some readers will be sufficiently intrigued by the connections that we have sought to establish 
in this paper to go on and grapple with them on their own. 
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